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Abstract 

One method of gaining some insight into the motion of particles 
in a medium with topological defects (e.g., electrons in a dislocated 
metal) is to look at the geodesies of the medium around the defect. 
In this work the Hamilton-Jacobi equation for the geodesies in a con- 
tinuous medium containing a torsional defect, an edge dislocation, 
is solved by using perturbation theory to first order in the Burgers 
vector. 
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Topological defects in solids have been described by Riemann-Cartan ge- 



have shown the equivalence between three-dimensional gravity with tor- 
sion and the theory of defects in solids. Based in this formalism, the non- 
relativistic quantum mechanics of electrons or holes around disclinations || |J 
and some properties of the quantum electromagnetic field also around discli- 
nations H H have been studied. Disclinations, defects associated with the 
breaking of rotational symmetry, are common in lower dimensional systems 
like liquid crystals |7j and an important ingredient in the making of ge- 
ometrically frustrated amorphous solids ||. They do not appear in bulk 
crystals,though, due to the very high elastic energy involved in their for- 
mation. Dislocations, associated with translational symmetry breaking, do 
appear in ordinary 3D crystals, but the asymmetry the! y provoke in the 
medium makes it difficult to study dynamical properties of objects moving 
in their midst. By presenting the geodesies around a dislocation, this work 
contributes to the visualization of the dynamics of such objects and hopefully, 
will provide a means of obtaining some insight into the motion of particles 
and the geometrical optics in dislocated media. 

We are concerned here with the geometrical structure of a crystalline 
medium that contains a single edge dislocation. An edge dislocation may be 
formed || by placing a pair of opposing parallel disclinations next to each 
other as a dipole (Fig. 1). Each disclination is made by either removing or 
inserting a wedge of material. In the continuum limit, at large distances from 
the defect, the medium is described by the metric [[J 



in cylindrical coordinates. The wedge angles are given by ±m and the discli- 
nations are at a distance h apart. The Burgers vector b = —mhx is therefore 
along the x-axis and the defect itself is along the z-axis. Since the defect 



ometry since the early 1950's 



[0]. More recently, Katanaev and Volovich || 




(1) 
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preserves translational symmetry along the z-axis, from now on we will be 
working in a z=constant surface. 

The geodesies may be obtained by requiring that the distance between 
two given points in the manifold be a minimum. That is, if 7 is a curve 
parameterized by t joining the points at t — to and at t — t\ in the manifold, 

ds 2 = gijdx 1 dx 1 (2) 

leads to 

sh] = hg^x^dt. (3) 

Jto 

Requiring s[~f] to be minimum, determines the path 7. The geodesic 7 is 
then a solution to the equation 

d 2 x % : dx j dx k 

where the Christoffel symbols of the second kind, T l - k) are 

V i 1 is ( 9 9sj dg sk _ dgjk\ 
lh 2 y \dx k dxi dx s ) K ' 

and gijgv = 1. 

In the language of Lagrangian mechanics, spy] is the action functional 
and the integrand in equation (3) is the Lagrangian. This procedure, when 
applied to the metric (1) leads to a cumbersome pair of coupled differential 
equations of second order. We follow instead the Hamilton- Jacobi method 
described in 0. But first we rewrite the z=const. section of metric (1) in 
Cartesian coordinates and, since we are interested in effects far away from 
the defect (i.e. r >> h), neglect the term in h 2 /r 2 : 



ds 2 =[l- h —^—\{dx 2 + dy 2 ). (6) 
7T x z + y l 1 



Justified by this same argument we will perform below perturbation theory 
in powers of b. 
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It is well known |TIJ that the Newtonian-like action 

W[l\ = \ t g^Vdt (7) 



also leads to the equations of geodesies. Using this action we proceed to 
write Hamilton- Jacobi equation 

dW TT (dW\ 

where H = L = ^gijX l & is the Hamiltonian and L is the Lagrangian. Now, 
writing H in terms of the moments Pi = §§r = f^" we get 

/ b y \dW (dW\ 2 ( dW\ 2 



It follows that 

W = c x t + W(x,y) (10) 

implying 

(dW\ 2 (dW\ 2 ( b y \ 

Next, we write the perturbation series for W 

W = W + bW l + b 2 W 2 + ... (12) 
Substituting (12) into (11) and equating terms of same order in b we find 

1 ] + — - +2 Cl = (13) 



and 



dx I \ dy 



dW dWt dW dW 1 = Cl y 

'a.. a.. .9. i v / 



dx dx dy dy ir x 2 + y 2 
up to first order corrections. 
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Equation (13) is easily solved by the ansatz 

W = c 2 x + c 3 y, 

where the constants ci, c 2 and c 3 are such that 

1 



^(4 + 4)- 



Now we turn to equation (14) which becomes 

C2— h c 3 



7T x 2 + y 2 



(15) 



(16) 



(17) 



<9x dy 

This equation is solved by the method of characteristics [IT]: first differenti 
ate W\ with respect to the auxiliary variable u, 



dW 1 (dx\ (dW x 



du 



_ + dy\ f dW 1 

duj\dx J \du J \ dy 



(18) 



then by comparing (17) and (18) make the identifications 



dWi ci 



dx 



dy 



du 



tc x 2 + y 2 ' du 2 ' du 



c 3 - 



(19) 
(20) 



It follows that 

x = c 2 u + ti, and y = c 3 u, 

where ti is another auxiliary variable. Now we substitute (20) in the first 
equation in (19) obtaining 



C1C3 



u 



(21) 



du 7T (c 2 U + ti) 2 + (c 3 «) 2 

Integration with respect to u of the above equation, followed by substitution 
of u = y/c 3 and t\ — x — —y (from (20)) and use of (16) yields 

C3 



W w at+Wo+bWx 



{c 2 +c 3 )t+c 2 x+c 3 y— - — m(x +y )H arctan 

47T ' 1_ 



2tt 



c 2 x + c 3 y 

c 3 x - c 2 y y 
(22) 



The corresponding solution (up to a constant A) of the equations of motion 
is obtained by differentiation with respect to c 2 and C3: 




- c 2 t + x + — arctan ( — ) + — ( 2 - 2 ) = (23) 

and 

-c 3 t + y-^\n(x 2 + y 2 )-^- [ ) = °- ( 24 ) 




Inserting the missing constant, eliminating t between the above equations 
and making the substitution B = c 3 /c 2 , we obtain 



y = A + Bx+^- 



- ln(x 2 + y 2 ) + B arctan ( - + — J + I 
2 \Bx-y 



(25) 



which gives the correct limit in the absence of the defect (6 = 0). A first 
iteration of this equation (still keeping to first order in b) allows us to plot 
the geodesies, which are shown in Figs. 2-4. Each figure shows a bundle of 
geodesies parallel to each other at infinity, approaching and being deflected 
by the dislocation located at the origin of the plots. 

Due to the particular choice of coordinates, plotting geodesies that asymp- 
totically are of the kind x=const. (i.e. vertical), is tricky. To do so, we first 
divide equation (25) by B and then take the limit B — > 00. We are left with 

C = x + — arctan (-) , (26) 



2tt \x 

where C is constant. As before, this gives the correct vertical geodesies in 
the limit b = 0. Inversion of the above equation gives 



y = x tan 



(27) 



Because of the periodicity of tan we need only to choose a convenient range 
for 6 = ^r(C — x) in order to plot the geodesies. Since our approach is only 
valid away from the origin, the range — 7r/2 < 6 < n/2 is excluded. We use 
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then the ranges — 3n/2 < 9 < —n/2 and n/2 < 9 < 3n/2. The result is 
shown in Fig. 5. 

In conclusion, the solution to the equation of geodesies in a medium 
containing one edge dislocation has been found by using perturbation theory 
to first order in the Burgers vector. A single dislocation produces a dilatation 
that breaks the translational symmetry of the crystalline lattice causing a 
deformation of the geodesies around the defect as it is clear from Figs. 2- 
5. This gives an intuitive picture to what happens to (classical) electrons 
moving in a dislocated metal or to the propagation of rays in a dislocated 
continuous medium, for example. 
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Figure Captions 

Figure 1: Disclination dipole making up an edge dislocation. 
Figure 2: Geodesies around edge dislocation, B — — 1. 
Figure 3: Geodesies around edge dislocation, B = 0. 
Figure 4: Geodesies around edge dislocation, B = 1. 
Figure 5: Geodesies around edge dislocation, B = oo. 
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Fig. 1 
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Fig. 4 




Fig. 5 




